Phase transitions in finite systems = topological peculiarities 
of the microcanonical entropy surface 
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It is discussed how phase transitions of first order (with 
phase separation and surface tension), continuous transitions 
and (multi)-critical points can be defined and classified for 
finite systems from the topology of the energy surface e'^'^'^' 
of the mechanical N-body phase space or more precisely of 
the curvature determinant D{E, N) = d"^ S/dE^ *d'^ S/dN^ - 
(d^S/dEdN)^ without taking the thermodynamic limit. The 
first calculation of the entire entropy surface S{E, N) for a 
q = 3-states Potts lattice gas on a 50 * 50 square lattice is 
shown. There are two lines, where S{E, N) has a maximum 
curvature ~ 0. One is the border between the regions in 
{E,N} with D{E,N) > and with D{E,N) < 0, the other 
line is critical starting as a valley in D{E, N) running from 
the continuous transition in the ordinary q = 3-Potts model, 
converting at Pm into a flat ridge/plateau (maximum) deep 
inside the convex intruder of S{E, N) which characterizes the 
first order liquid-gas transition. The multi-critical point Pm 
is their crossing. 
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Boltzmann's gravestone has the famous epigraph: 



k * InW 



which puts thermodynamics on the ground of mechan- 
ics [0,1). It relates the entropy S to the volume 
W{E, N, V) ^ Se* trS{E - Hn) of the energy (E) sur- 
face of the N-body phase space at given volume {V), the 
microcanonical partition sum. Here de is a suitable small 
energy constant, Hn is the A^-particle Hamiltonian, and 
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The set of points on this surface defines the microcanon- 
ical ensemble [ME). 

Today conventional thermodynamics is based on the 
canonical statistical mechanics as introduced by Gibbs 
In the thermodynamic limit ThL (iV— >oo|jv/y=const) 
the canonical ensemble {CE) is equivalent to the funda- 
mental {ME) if the system is in a pure phase and the 
ThL exists. 

The fundamental difference of microcanonical ther- 
modynamics (MT) to conventional thermodynamics is 
that no non-mechanical quantities like temperature, heat, 
pressure have to be introduced a priori. 

The link between ME and CE is established by Laplace 
transform. E.g. the usual grand canonical partition sum 
is the double Laplace transform of ME: 



Z{T,ijL,V)= fr dE dN e-^^'f'^^/^trSiE - Hn). (2) 



This excludes all inhomogeneous situations, especially 
phase separations. There the entropy is non-extensive 
and the CE contains several Gibbs states at the same 
temperature. Consequently, the statistical fluctuations 
do not disappear in the CE even in the thermodynamic 
limit. This is the reason why Gibbs himself excluded 
phase separations in chapter VII of |^ in a footnote, page 
75. At phase transitions the ME and the CE describe dif- 
ferent physical situations. If one combines a small system 
of water at the specific energy ei of boiling water with 
a large heat bath at 100°C it will remain at 100°C but 
may convert into steam in 50% of the cases. I.e. the fun- 
damental assumption used e.g. by Einstein |^ that our 
"system changes only by infinitely little" does not hold. 

It is important to notice that Boltzmann's and also 
Einstein's formulation allows for defining the en- 

tropy by Smicro ■— ln[W{E, N, V)] (in the following we 
use Smicro for S{E, N) if it is not clear) as a single valued, 
non-singular, in the classical case differentiable, function 
of all "extensive", conserved dynamical variables. No 
thermodynamic limit must be invoked and the theory ap- 
plies to non-extensive systems as well. Of course this is 
achieved by avoiding Gibbs-states, "equilibrium states" 
or "most random" states [Q. On the other hand fluc- 
tuations become then important and must be simulated 
by Monte Carlo methods. The microcanonical ensemble 
is the entire microcanonical N-body phase space without 
any exception. In MT the entropy is not "a measure of 
randomness" [Q, it is simply the volume e^^^'^'^^ of the 
energy surface. The latter point is extremely important 
as it allows to address even thermodynamically unstable 
systems like collapsing gravitating systems (for a recent 
application of MT to thermodynamical unstable, collaps- 
ing systems under high angular momentum see |^ ) . In so 
far it is the most fundamental formulation of equilibrium 
statistics From here the whole thermostatics may 

be deduced. MT describes how e'^^-^'-'^'^) depends on the 
dynamically conserved energy, number of particles etc.. 
Of course we must assume that the system can be found 
in each phase-space cell of e"^ with the same probability. 

Following Lee and Yang [|| phase transitions are in- 
dicated by singularities in Z{T, fi,V). Singularities of 
Z{T, n,V), however, can occur in formula (^) in the ther- 
modynamic limit only {V —^ oo\N/v=g,E/N=e)- Foi' finite 
volume ^(T, /i, V) is a finite sum of exponentials and ev- 
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erywhere analytical. Only at points where S{E,N) has 
a curvature > will the integral eq.|^ diverge in the ther- 
modynamic limit. In these points, the Laplace integral]^ 
does not have a stable saddle point. Here van Hove's con- 
cavity condition ||] for the entropy S{E, N, V) of a stable 
phase is violated. Consequently we define phase transi- 
tions also for finite systems topologically by the points of 
non-negative curvature of the entropy surface S{E, N, V) 
as a function of the mechanical, conserved "extensive" 
quantities like energy, mass, angular momentum etc.. 

Experimentally one identifies phase transitions of 
course not by the singularity of Z{T, /x) but by the inter- 
faces separating coexisting phases, e.g. liquid and gas, 
i.e. by the inhomogeneities. The interfaces have three 
effects on the entropy: 

1. There is an entropic gain by putting a part (A^i) 
of the system from the majority phase (e.g. hquid) 
into the minority phase (bubbles, e.g. gas), but this 
is connected with an energy-loss due to the higher 
specific energy of the "gas" -phase, 

2. an entropic loss proportional to the interface 
area by the correlations between the particles in 
the interface, leading to the convex intruder in 
S{E, N, V) and is the origin of surface tension pO[ , 

3. and an additional mixing entropy for distributing 
the A^i-particles in various ways into bubbles. 

At a (multi-) critical point two (or more) phases become 
indistinguishable and the interface entropy (surface ten- 
sion) disappears. 

Microcanonical thermostatics was introduced in great 
detail for atomic clusters and nuclei in ||l^ . In two further 
papers we showed how the surface tension can quantita- 
tively be determined from the microcanonical S{E, N, V) 
for realistic systems like small liquid metal drops ||l2|,|l3| . 
Needless to say that only by our extension of the defini- 
tion of phase transitions to finite systems it make sense 
to ask the question: "How many particles are needed to 
establish a phase transition?" The results we are going 
to show here will again demonstrate that often N does 
not need to be large in order to see realistic values (close 
to the known bulk values) for the characteristic parame- 
ters. Other examples where shown in our earlier papers 

In the following we investigate the 3-state Potts lattice- 
gas model on a 2-dim = 50^ square lattice. We will 
see how the total microcanonical entropy surface S{E, N) 
decovers even the most sophisticated thermostatic fea- 
tures as first order phase transitions, continuous phase 
transitions, critical and multicritical points even for fi- 
nite systems and non-extensive systems. This demon- 
strates that the microcanonical statistics is in contrast 
to Schrodinger's claim Q quite able to handle not only 
gases but also phase transitions and critical phenomena. 



Important details like the separation of phases and the 
origin of surface tension can be treated. The singularity 
of the canonical partition sum at a transition of first or- 
der can be traced back to the loss of entropy due to the 
correlations between the surface atoms at phase bound- 
aries pl| . 




FIG. 1. Entropy s = -^S',„>cro(e=£'/L2,n=iV/L2). The 
grey levels are as in fig.^ white regions: concave, pure phases 
(in APmC ordered, "solid", in CPmB disordered "gas"); 
black: convex, phase separation ("liquid-gas"); and light grey 
strips: critical branches; crossing: multicritical point Pm- If 
one plots Smicro VS. /3 and like in CT the two wings PmA 
and PmB are mapped onto one-another in the figure H below. 
The black regions of the intruder get folded in between (see 
figi- 

Briefly, a few words about our method which will be 
published in detail in |^: The Hamiltonian is H = 
~^J2ij^<Ti,<yj nearest neighbors). We covered all 
space {E = e*L^,N = n*L^} by a mesh with about 1000 
knots with distances of Ae = 0.04 and An — 0.02. Due to 
our limited computational resources (DEC- Alpha work- 
station) we could not use a significantly denser mesh. 
At each knot {ei,nk} we performed by microcanonical 
simulations (« 2 * 10^ events) a histogram for the prob- 
abilities P(ei,n}S) for the system to be in the narrow re- 
gion [Ei ± 4) * {N}^ ± 4) of phase space. Local deriva- 
tives (3 = {dS{E,N)/dE)^, p^i ^ - {dS{E, N)/dN)j^ 
in each histogram give the "intensive" quantities, so that 
the entire surfaces of S{E,N), I3{E,N), (3^i{E,N) can 
be interpolated. The first derivatives of the interpolated 
(smoothed) (3{E, N) and f}fi{E, N) give the curvatures. 

The figure ^ shows some of our recent results 
for S{E/L^,N/L^) for the case of the diluted 
(7 = 3 Potts model. Grid lines are in direc- 
tion [E - EQ{N)]/[E„,ax{N) ~ Eq{N)] =const. resp. 
N/E^ =const.. The black region is the intruder at the 
first-order condensation transition ("liquid-gas coexis- 
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tence") with positive largest curvature of S{E,N). This 
corresponds to the similar region in the Ising lattice gas, 
respectively the original Ising model as function of the 
magnetization. At the light grey strip S{E, N) is critical 
with vanishing largest curvature. The line from point C 
over the multicritical point P„i to D corresponds from 
C to Pm to the familiar continuous transition in the or- 
dinary g = 3 Potts model. At Pm this line crosses the 
rim of the intruder from A to B which is the border of 
the first order transition. This crossing determines the 
multicritical point Pm quite well at (3m = 1-48 ± 0.03, 
Pf^m = —2.67 ± 0.02 or Cm —1, ~ 0.7. From here 
the largest curvature starts to become > 0. Naturally, Pm 
spans a much broader region in {e, n} than in {/?, 
remember here S{E, N) is flat. 



D{E,N) = 




FIG. 2. Contour plot of the entropy SmicroiP, PfJ,). One 
clearly sees the critical line where Smicro{l3, PfJ,) becomes 
tripelvalued and/or at the ordinary Potts continuous transi- 
tions dS/dp ~ ~oo. Here several points of S{E, N) from the 
entropy surface are folded onto one-another. That is the ori- 
gin of the thick multi- valued ("critical") line in Smicro{P, P/J-)- 

If one plots the entropy Smicro{P, /?^) as function of the 
"intensive" variables |3^i = -dS/dN and (3 = dS/dE, 
we obtain picture |[ This corresponds to the conven- 
tional grand-canonical representation if we would have 
calculated the grand canonical entropy from the Laplace 
transform Z{T,ijl,V), eq.^. As there are several points 
Ei,Nk with identical (3,(3^, Smicro{(3, (3^) is a multival- 
ued function of (3,(3^. Here the entropy surface S{E,N) 
is folded onto itself see fig.^ and in fig.|| these points show 
up as a black critical line (dense region) . The backfolded 
branches of S{E, N) are jumped over in eq. || and get 
consequently lost in Z{T,^). This demonstrates the far 
more detailed insight one obtains into phase transitions 
and critical phenomena by microcanonical thermostatics 
which is not accessible by the canonical treatment. 

In figure the determinant of curvatures of S{E, N): 



dW dNdE 
OEdN dW 



(3) 



is shown. On the diagonal we have the ground-state of 
the 2-dim Potts lattice-gas with e = —2n, the upper-right 
end is the complete random configuration (not shown), 

9 2 

with the maximum allowed excitation erand — In 
the upper right (white) D > 0, both curvatures are nega- 
tive. In this region the Laplace integral eq.^ has a stable 
saddle point. This region corresponds to pure phases. 



beta 




-4.25 '^-^^ 



-3.25 



-2.75 



bemu 



FIG. 3. SmicroiP, PfJ-), the same as figj^. The grey lines are 
for (3 =const. The backbending of the intruder left of Pm 
which is hidden behind the thick critical line in fig.^ is clearly 
seen. This corresponds to the whole intruder of phase separa- 
tion (black and grey in fig.^). The lines connect the points of 
ds/d{l3^)\p — — oo and thus correspond to the critical lines in 
fig.^ The positions of A,B,D,C are only roughly indicated. 

In the light gray region we have D ^ 0. This is the 
critical region. Here the largest eigenvalue of D is 0. Two 
branches cross here: One goes « parallel to the ground 
state {E oc -27V) from A to B. This is a rim in D{E, N), 
the border line between the region with D{E, N) > 0, 
and the region with D{E, N) < (black) where we have 
the first order liquid — gas transition of the lattice-gas. 
The Laplace integral (||) has no stable saddle point and in 
the ThL the grand canonical partition sum (|^) diverges. 
Here we have a separation into coexisting phases, e.g. 
liquid and gas. Due to the surface tension or the negative 
surface entropy of the phase boundaries, S{E,N) has a 
convex intruder with positive largest curvature. 

The other branch from C to Pm is a valley in D{E, N). 
Here the largest curvature of S{E, N) has a local mini- 
mum and D ^ (it would be D = with a higher preci- 
sion of the simulation), running from the point (near C) 
of the continuous phase transition at n = 1 and e = —1.57 
of the ordinary q = 3-Potts model downwards to Pm- It 
converts below the crossing point Pm into a flat ridge 
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inside the convex intruder of the first order lattice-gas 
transition. The area of the crossing of the two critical 
branches CPmD and AP^B is the multi- critical region 
Pm of the q — 3 Potts lattice gas model. 




-1.75 -1.5 -1.25 -1 -0.75 -0.5 -0.25 
epsilon 

FIG. 4. Contour plot of the determinant of curvatures 
D{€—E/L'^,n=N/L^). In the two light gray strips we have 
D{e,n) ~ 0. Here the transition is critical (continuous). The 
crossing point Pm indicates the multicritical region. The two 
inlets show cuts along the line /<->(/ (e — — 2n-|-0.28) through 
the intruder with first order transition and negative D (pos- 
itive largest curvature). To understand the dramatic loss 
of information in the CT: The whole black and grey region 
(about half of the entire phase space !) gets lost, see also fig^ 
The other along m<^n near to the ordinary q=3 Potts model 
(n = 0.88) shows to the left the ordered phase with positive 
D, then the deep valley with nearly vanishing curvature and 
to the right the disordered ("gas") phase. 

Conclusion: Microcanonical thermostatics (MT) de- 
scribes how the entropy S{E, N) as defined entirely in 
mechanical terms by Boltzmann depends on the con- 
served "extensive" mechanical variables: energy E, par- 
ticle number TV, angular momentum L etc. This allows to 
study phase transitions also in small and in non-extensive 
systems. If we define phase transitions in finite systems 
by the topological properties of the determinant of cur- 
vatures D{E,N) (eq.||) of the microcanonical entropy- 
surface S{E, N): a single stable phase by D{E, A^) > 0, a 
transition of first order with phase separation and surface 
tension by D{E,N) < , a continuous ("second order") 
transition with D{E, N) = 0, and a multi- critical point 
where more than two phases become indistinguishable 
by the branching of several lines with D{E, N) — 0, then 
there are remarkable similarities with the corresponding 
properties of the bulk transitions. 

The advantage of MT compared to CT is clearly 
demonstrated: About half of the whole phase space, the 



intruder of S{E, N) or the non-white region in fig.^, gets 
lost in conventional canonical thermodynamics. Without 
any doubts this contains the most sophisticated physics 
of this system. Due to limited computer resources this 
could be demonstrated with only limited precision. We 
are convinced our conclusions will be verified by more 
extensive - and more expensive - calculations. 
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